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Abstract 

In this paper, we introduce the generalized Gompertz-power series class of distribu¬ 
tions which is obtained by compounding generalized Gompertz and power series distribu- 
tions. This comp oundi ng procedure follows same way that was previously carried out by 


Silva et al 


(2013) and 


Barreto-Souza et al. 


(|201lll in introducing the compound class of 


extended Weibull-power series distribution and the Weibull-geometric distribution, respec¬ 
tively. This distribution contains several lifetime models such as generalized Gompertz, gen¬ 
eralized Gompertz-geometric, generalized Gompertz-poisson, generalized Gompertz-binomial 
distribution, and generalized Gompertz-logarithmic distribution as special cases. The haz¬ 
ard rate function of the new class of distributions can be increasing, decreasing and bathtub¬ 
shaped. We obtain several properties of this distribution such as its probability density 
function, Shannon entropy, its mean residual life and failure rate functions, quantiles and 
moments. The maximum likelihood estimation procedure via a EM-algorithm is presented, 
and sub-models of the distribution are studied in details. 

Keywords: EM algorithm, Generalized Gompertz distribution, Maximum likelihood es¬ 
timation, Power series distributions. 

1 Introduction 


The exponential dist ribution i s com monly used in many applied problems, particularly in life¬ 


time data analysis (jLawlessl . 


20031 ). A generalization of this distribution is the Gompertz 


distribution. It is a lifetime distribution and is often applied to describe the distribution of 
adult life spans by actuaries and demographers. The Gompertz distribution is considered for 
the analysis of sur vival in some sciences such as biology, gerontology, computer, and marketing 
science. Recently, iGuota and Kundul (|l999l i defined the generalized exponential distribution 
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and in similar manner. lEl-Goharv et alJ (|2013l i introduced the generalized Gompertz (GG) dis¬ 
tribution. A random variable X is said to have a GG distribution denoted by GG(a,/ 3 , 7 ), if 
its cumulative distribution function (cdf) is 


G(x} = [1 — a,/3 > 0, 7 > 0; x > 0. 

and the probability density function (pdf) is 


( 1 . 1 ) 


( 1 . 2 ) 


The GG distribution is a flexible distribution that can be skewed to the right and to the 
left, and the well-known distributions are special cases of this distribution; the generalized 
exponential proposed by 


Gupta and Kundul (jl9^ 1 when 7 —)• O'*', the Gompertz distribution 


when 0 = 1, and the exponential distribution when a = 1 and 7 —?■ O'*'. 

In this paper, we compound the generalized Gompertz and power series distributions, 
and introduce a new class of distribution. This procedure follows similar way that was previ- 
ously carried out by some aut hors: The exponential-power series distribution is introduced by 


Ghah 


^nnJ 


candi and Ganialil (1200911 w 


Adamidis and Loukas 


(Tahmasbi and Rezaei. 


by 


2008 


rich is concluded the e xponential 


1998|), exponential-Poisson (|Kua . 120071 1. and exponential-logarithmic 


-geometric (jAdamidis et al. 


distri 


rutions; the Weibull-power series distributions is introduced 


Morais and Barreto-Souzal (j201lh and is a generalization of the exponential- power series dis¬ 


tribut ion; the generalized exponential-power series di stribution is in t roduced by 


2012l i which is concluded the Poisson-expone ntial (iGancho et al 


201 ih complementary ex ponential-geometric (jLouzada et al 


exponentia’ 


-power series (IFlores et al 


tributions (jMahmoudi and .Tafari 


201 4l i. 


2011 


Mahmoudi and Jafari 


Louzada-Neto et al. 


201 ih . and the complementary 


2013 ) distributions; linear failure rate-power series dis- 


The remainder of our paper is organized as follows: In Section [2l we give the probability 
density and failure rate functions of the new distribution. Some properties such as quantiles, 
moments, order statistics. Shannon entropy and mean residual life are given in Section O In 
Section [H we consider four special cases of this new distribution. We discuss estimation by 
maximum likelihood and provide an expression for Fisher’s information matrix in Section [5j A 
simulation study is performed in Section [H An application is given in the Section [71 










































































2 The generalized Gompertz-power series model 


A discrete random variable, is a member of power series distributions (truncated at zero) if 
its probability mass function is given by 


Pn = P{N = n) = 


anO^ 

c{e) 


, n = 1,2 ,..., 


( 2 . 1 ) 


where a„ > 0 depends only on n, C{9) = and 6 E (0, s) (s can be oo) is such that 

C{9) is finite. Table [T] summarizes some particular cases of the truncated (at zero) power series 


distributions (geometric, Poisson, 
series distribution can be found in 


ogarit 


i mic a nd binomial). Detailed properties of power 


Noack ( 1950 b Here, C'{9), C''{9) and C'''{9) denote the 


first, second and third derivatives of C{9) with respect to 9, respectively. 


Table 1: Useful quantities for some power series distributions. 


Distribution 

an 0(9) C'(9) C"(9) 0'''(9) s 

Geometric 

Poisson 

Logarithmic 

Binomial 

1 9(1-9)-^ (1 - 6»)-2 2(1 - 6»)-3 6(1-6 »)-^ 1 

n\~^ — 1 e® oo 

n-^ -log(l-0) (1-9)-^ (1-0)”^ 2(1-9)-^ 1 

('in\ 1 a\m 1 m m(m-l) m(m-l)(fc-2) 

\nJ G “T^ 


We define generalized Gompertz-Power Series (GGPS) class of distributions denoted as 
GGPS(a, /3, 7 , 9) with cdf 


T(x) f^ ^r.{eG{x)r C{9G{x)) C{9n ^ ^ ^ 


n=l 


C{9) 


C{9) C{9) 


( 2 . 2 ) 


where t = I — e t" 


. The pdf of this distribution is given by 


(2.3) 


This class of distribution is obtained by compounding the Gompertz distribution and power 
series class of distributions as follows. Let be a random variable denoting the number of 
failure causes which it is a member of power series distributions (truncated at zero). For 
given N, let Xi,X 2 , ■ ■ ■ ,Xn' be a independent random sample of size N from a GG{a, 13,'y) 
distribution. Let = maxi<j<jv A"*. Then, the conditional cdf of | A = n is given by 

_ [1 _ p-f 


which has GG(na,/3,7) distribution. Hence, we obtain 


P(A(7 v) < x,N = n) = 


ani9G{x))^ ^ an9^ _ £(e7^_i) „ 

0(9) Ci9)'- ‘ ■ 


















Therefore, the marginal cdf of has GGPS distribution. This class of distributions can 

be applied to reliability problems. Therefore, some of its properties are investigated in the 
following. 


Proposition 1. The pdf of GGPS class can be expressed as infinite linear combination of pdf 
of order distribution, i.e. it can he written as 

OO 

f{x) = '^Png{n){x;na,P,-f), (2.4) 

71=1 

where g(„)(x; no,/3, 7 ) is the pdf of GG{na, 

— §. — V) 

Proof. Gonsider t = 1 — e t" ^ h So 

^ /7 

= = '^Png{n){x;na,ld,j). 

71 = 1 ^ 71=1 

□ 


Proposition 2. The limiting distribution 0 / GGPS(a,/3, 7 , 0) when 0 —)• O'*" is 


lim F{x) = [1 — e i)jca 

6 ^- 0 + 


which is a GG distribution with parameters ca, fd, and 7 , where c = min{n G N : o„ > 0}. 


Proof. Gonsider t = 1 — e 




. So 


lim F(x) 
9 ^ 0 + 


C(At"' 

el^+~C{9) 


lim 

A->-0+ 


OO 


71=1 


OO 


E anO^ 

71=1 


= lim 

0 ^. 0 + 


OO 

aci“+ E 

n=c+l 

OO 

flc + E CLnO"-~'^ 

n=c+l 


= F 


□ 


Proposition 3. The limiting distribution 0 / GGPS(a, /?, 7 , 0) when 7 —?■ O"*" is 


lim F(x) 

7 ^ 0 + 


C( 0(1 - 

W) 


i.e. the cdf of the g e neral ized exponential-power series class of distribution introduced by 
Mahmoudi and .Jafari f20lA ). 












Proof. When 7 ^ 0 ^, the generalized Gompertz distribution becomes to generalized exponen¬ 
tial distribution. Therefore, proof is obvious. □ 


Proposition 4. The hazard rate function of the GGPS class of distributions is 


h{x) 


eal3e'^^{l-t)t^-^C'{et°^) 

C{e) - 


(2.5) 


ILi. 

where t = 1 — e r 


(eT'^-l) 


Proof. Using ( 12 . 21 ) . (|2.3I) and definition of hazard rate function as h{x) = f{x)/{l — F{x)), the 
proof is obvious. □ 


Proposition 5. For the pdf in (12.3p . we have 


lim /(x) 
X—> 0 + 


C (0)9(3 
C(9) 


0 < a < 1 
0 = 1 
a > 1, 


lim f{x) = 0 . 

X^OO 


Proof. The proof is a forward calculation using the following limits 


lim 

a:->- 0 + 


'00 0 < a < 1 

< 1 0 = 1 
_ 0 o > 1 , 


lim F = 0 , lim t = 1 . 
a;->- 0 + x^oo 


□ 


Proposition 6. For the hazard rate function in (|2.5p . we have 


lim h{x) 

3;^0+ 


00 

c'(o)ep 

0 ( 9 ) 


0 < o < 1 
0 = 1 
o > 1 , 


lim h{x) 

x—>-oo 


00 7 > 0 
/3 7 —)• 0 


Proof. Since lim 3 ,^o+(l “ F{x)) = 1, we have lim 2 ,_^Q+ h{x) = lim 2 ,_^Q+ f{x). 
For lima;_>.oo h{x), the proof is satished using the limits 


lim c'{er) = c'(e), lim = i, 

X—)-oo x^oo 

lim = lim ^ f 7 > 0 

x^co C(9) — C(Otr x^oo 0l3aC'{6)e^^{l — t) y 9aC'(9) 7 0. 


□ 

As a example, we consider C {9) = 6 + 9“^^. The plots of pdf and hazard rate function of 
GGPS for parameters /3 = 1,7 = .01,0 = 1.0, and o = 0.1,0.5,1.0, 2.0 are given in Figured] 


This pdf is bimodal when o = 2.0, and the values of modes are 0.7 and 3.51. 
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Figure 1: Plots of pdf and hazard rate functions of GGPS with C {9) = 6 -\- 6“^^. 


3 Statistical properties 

In this section, some properties of GGPS distribution such as quantiles, moments, order statis¬ 
tics, Shannon entropy and mean residual life are obtained. 


3.1 Quantiles and Moments 


The quantile q of GGPS is given by 


X 




^ ^ 9 ^ 


0 < O' < 1, 


where G~^{y) = ^log[l — ) ] and C'“^(.) is the inverse function of C{.). This result 

helps in simulating data from the GGPS distribution with generating uniform distribution data. 

For checking the consistency of the simulating data set form GGPS distribution, the his¬ 
togram for a generated data set with size 100 and the exact pdf of GGPS with C (9) = 9 + 9'^^, 
and parameters a = 2, /3 = 1, 7 = 0.01, 9 = 1.0, are displayed in Figure [2] (left). Also, the 
empirical cdf and the exact cdf are given in Figure [2] (right). 

Consider X ~ GGPS(a,/3, 7 ,0). Then the Laplace transform of the GGPS class can be 
expressed as 


L{s) = E{e-^^) = Y^P{N = n)Ln{s), 

n=l 


where Ln{s) is the Laplace transform of GG(na,/ 3 , 7 ) distribution given as 

r+oo 


LJs) = 


f 




= na/J 


I 




E 

j=0 




(3.1) 
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Figure 2: The histogram of a generated data set with the pdf (left) and the empirical cdf with 
cdf (right) of GGPS distribution. 


na — 1 

j 


{-lye 


na/3j; ("Y‘) / 


{'1-s)x^=y(j+l)e^- 


dx 


T{k + l) 


nafd^Yl 


na 


-l\ (-l)^'+^ef(^'+')[^(j + l)] 


J 


T{k + l)(s — 7 — 'jk) 


s > 7- 


(3.2) 


j=0 k=0 

Now, we obtain the moment generating function of GGPS. 

OO 

Mx{t) = E{y^) = Y,P{N = n)Ln{-t) 


n=l 


n=l ' ' k=0j=0 


- n(’^7')(-l)^+"+7^^^+'7f(j + 1))* 


T{k + l)(t + 7 + 'jk) 

S.I 


- “ iV(^“-i)(-l)Wie7(^'+i)(^(j + 1))^ 

fc=0 j=0 


(3.3) 


r(A: + l)(t + 7 + 7A;) 

where is a random variable from the power series family with the probability mass function 
in (EU) and E]y[U] is expectation of U with respect to random variable N. 

We can use Mx{t) to obtain the non-central moments, = ElX'"]. But from the direct 
calculation, we have 


n=l ' ' k=0 j=0 


al3Ex[^'^ 

k=0 j=0 


r(k -|- 1)(7 + jky~^^ 

^(7Va-l)(_^)j+fc+r+lgf (i+l)(^(^- ^ + 1) 

r(k + i)(j + jky+^ ^ 


(3.4) 





























Proposition 7. For non-central moment function in\3.4[ we have 


lim fij. = E[Y^], 
e^o+ 


where Y has GG(ca, (3, 7 ) and c = min{n G N : > 0}. 


Proof. If y has GG(ca,/ 3 , 7 ), then 


^ ^ hh r{k + i){^ + jky+^ 


Therefore, 


E ane^E[Yy 

lim fir = lim ^ - 

0 ^ 0 + 0 ^ 0 + ^ 

n=l 


= lim 
e^-o+ 


00 

a,E[Yy+ E 

n=cH-l 

00 

«c + I] anO^-^ 
n=c+l 


= E[Yy. 


□ 


3.2 Order statistic 


Let Xi,X 2 , ..., Xn be an independent random sample of size n from GGPS(a, (3, 7 , 9). Then, 
the pdf of the ith order statistic, say Xi-n, is given by 

f. (x) = _ — _ V ~Ml \ n-i 

h:n[ ) ^i_i)i^n-i)r^c{e) ^ ^ c{9) ^ ’ 

_ _§.^Qyx_-^\ 

where / is the pdf given in ( 12 .31) and t = 1 — e r 7 Also, the cdf of Xi^n is given by 


F ui G (-i)n 7 ‘) , c(t“) 

' (i-l)!(n-i)!^ k + i+l 7(9) 


k-\-i 


An analytical expression for rth non-central moment of order statistics Xi-n is obtained as 

/t _ 1 \ /r)\ r+°° 

fc-n-l-j-l / ^ ^\ ''' 


k^—Tt — 1~\~ 1 


= r 


E 


(_l)fc-n+*-l _ 1 


k=n—i+l 


[C{9yy \n — ij \k 


:)r 


x^-yc{d) - c{ee)fdx 


where S{x) = 1 — E{x) is the survival function of GGPS distribution. 












3.3 Shannon entropy and mean residual life 


If X is a no ne-negativ e cont inuous random variable with pdf /, then Shannon’s entropy of X 


is defined by 


Sha,Tinonl ( 19481 ) as 


^- 1-00 

H{f) = E[- log f{X)] = - f{x) log{f{x))dx, 

Jo 

and this is usually referred to as the continuous entropy (or differential entropy). An explicit 
expression of Shannon entropy for GGPS distribution is obtained as 


Hif) = E{- log[^e^^(e-f- e'f(^(1 - g-f 

— (a — l)ili[log(l — — ili[log(C''^0(1 — 

= - “ 71^1 + -Mxil) - - - {a - 1)Y^ P{N = n) f log{t)dt 

E{d) 7 7 ^ Jo 

oo „! 

— ^2 P{X = n) nvP'~^ \og{C'{0u))du 

n=i Jo 

= -l<>gl^|-Wi + -Afe(7)-- + 7^EN|7]-£„[A(Ar,())], (3,5) 

0 ( 0 ) 7 7 a A 

where A{N, 0) = Nu^~^ \og{C'{6u))du, A is a random variable from the power series family 

with the probability mass function in (12.ip . and Ex\U] is expectation of U with respect to 
random variable N. In reliability theory and survival analysis, X usually denotes a duration 
such as the lifetime. The residual lifetime of the system when it is still operating at time s, is 
Xs = X — s\X>s which has pdf 


f{x]s) 


fix) eg{x)C'{eGix)) 
l-F(s) C( 0 )-C( 0 G(s))’ - 


Also, the mean residual lifetime of Xg is given by 


mis) = E[X - s\X > s] 


C°°ix - s)fix)dx 
1-Fis) 
xfix)dx 
l-Eis) "" 
Cie)EN[Zis,N)] 
C(0)-C(0[l-e“^^"^"“^^]“) 


where Z(s,n) = xg^n)(x','na,/3,j)dx, and g(„)(x; na,/1, 7 ) is the pdf of GG(na,/3, 7 ). 















4 Special cases of GGPS distribution 

In this Section, we consider four special cases of the GGPS distribution. To simplify, we 
consider t = 1 — x > 0, and Aj = 


4.1 Generalized Gompertz-geometric distribution 


The geometric distribution (truncated at zero) is a special case of power series distributions 
with On = 1 and C{0) = (0 < 9 < 1). The pdf and hazard rate function of generalized 

Gompertz-geometric (GGG) distribution is given respectively by 

(l-0)a/3eT'"=(l-f)r-i 


Remark 4.1. Consider 


fix) = 

h{x) = 

fuix) = 


{9t^ - 1)2 

(l-0)a/3eT'^(l-t)t“-i 
(1 - 6 »t")(l - t“) ' 


X > 0, 

X > 0. 


X > 0, 


(4.1) 

(4.2) 

(4.3) 


((1 - 6 »*)t“ - 1)2 ’ 

where 9* = 1 — 9. Then fuix) is pdf for all 9* > 0 (see 
that when a = 1 and 7 —0+, the pdf of extended exponential geometric (EEC) distribution 


Marshall and Olkin. 


1991). Note 


l Adamidis et al\. 12003) is concluded from (14.311 . The EEC hazard function is monotonically 


increasing for 9* > 1; decreasing for 0 < 0* < 1 and constant for 9* = 1. 

Remark 4.2. If a = 9* = 1, then the pdf in (14.3p becomes the pdf of Gompertz distribution. 
Note that the hazard rate function of Gompertz distribution is h{x) = jde"^^ which is increasing. 

The plots of pdf and hazard rate function of GGG for different values of a, /3, 7 and 9* 
are given in Figure [3l 

Theorem 4.1. Consider the GGG hazard function in (SSI). Then, for a > 1, the hazard 
function is increasing and for 0 < a < 1 , is decreasing and bathtub shaped. 


Proof. See Appendix A.l. 

The first and second non-central moments of GGG are given by 


□ 


Em - - 0 E n."- E E +1))‘ 

k=0j=0 


n=l 


T{k + 1)(7 -I- 7A:)2 


E{X^) = 2a/3{l -9)Y^ n9^-^ EE 

k=0j=0 


n=l 


7L.(_l)i+fe+3efO+i)(|(j + l))fc 

r(A:-I- 1)(7 + 7 ^)^ 
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Figure 3: Plots of pdf and hazard rate function of GGG for different values a, /?, 7 and 9*. 


4.2 Generalized Gompertz-Poisson distribution 

The Poisson distribution (truncated at zero) is a special case of power series distributions 
with On = ^ and C{9) = — 1 {9 > Q). The pdf and hazard rate function of generalized 

Gompertz-Poisson (GGP) distribution are given respectively by 


f{x) 

h{x) 


9al3e^^-\l-t)e-^e^^\ x>0 

0a/3eT'*(l 

g6» _ g6>t“ ’ X > U. 


(4.4) 

(4.5) 


Theorem 4.2. Consider the GGP hazard function in ra- Then, for a > 1, the hazard 
function is increasing and for 0 < a < 1, is decreasing and bathtub shaped. 


Proof. See Appendix A. 2 . 


□ 


The hrst and second non-central moments of GGP can be computed as 




- - 7l,-(-l)^'+^ef(^'+^)(f(j + 1 ))'= 

^ (n- 1 )! ^ ^ 

=1 ^ ^ k= 0 j =0 


T{k -I- 1)(7 -I- 




r 7, .7 (i +1))* 

- 1 (™- 1 )! P{k + l){-i+ -ikf 
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Figure 4: Plots of pdf and hazard rate function of GGP for different values a, /3, 7 and 6. 

The plots of pdf and hazard rate function of GGP for different values of a, /3, 7 and 9 are given 
in Figured! 


4.3 Generalized Gompertz-binomial distribution 

The binomial distribution (truncated at zero) is a special case of power series distributions with 
an = (™) and C{6) = (0 + 1)™ — 1 (0 > 0), where m {n < m) is the number of replicas. The 
pdf and hazard rate function of generalized Gompertz-binomial (GGB) distribution are given 
respectively by 


fix) 

(g+oi _j_ 1 \m—1 

= x > 0 , 

(4.6) 

h{x) 

m0a/3eT'^(l - + 1 )”^-^ 

“ {9 + 1 )™ - ( 0 t" -Fl)™ , x> . 

(4.7) 


The plots of pdf and hazard rate function of GGB for m = 4, and different values of a, 
(3, 7 and 0 are given in Figure O We can find that the GGP distribution can be obtained as 
limiting of GGB distribution if mO —?■ A > 0, when m —?• 00 . 


Theorem 4.3. 


Consider the GGB hazard function in 


Then, for a > 1, the hazard 


function is increasing and for 0 < a < 1, is decreasing and bathtub shaped. 
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Figure 5: Plots of pdf and hazard rate function of GGB for m = 5, and different values a, (3, 
7 and 9. 


Proof. The proof is omitted, since 9 > 0 and therefore the proof is similar to the proof of 
Theorem 14.21 □ 


The first and second non-central moments of GGB are given by 


E(X) 


E(X^) 


fm] + 1 ))" 

(0 + 1 )”*-!^^ r(k + i)(j + jkf 

2a/3 ^ fm\ ^ ^ + 1))" 

^\nj r(k + l)(j + jk)3 


4.4 Generalized Gompertz-logarithmic distribution 


The logarithmic distribution (truncated at zero) is also a special case of power series distribu¬ 
tions with On = ^ and C(9) = — log(l — 9) (0 < 9 < 1). The pdf and hazard rate function of 
generalized Gompertz-logarithmic (GGL) distribution are given respectively by 


f(x) 

h{x) 


9ape^^{l-t)t^-^ 

— 1) log(l — 9) ’ 

9aPe^^{l-t)t^-^ 
{9t^ - l)log(^^)’ 


X > 0, 


X > 0. 


(4.8) 

(4.9) 
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Figure 6 : Plots of pdf and hazard rate function of GGL for different values a, /?, 7 and 9. 


The plots of pdf and hazard rate function of GGL for different values of a, j3, 7 and 9 are 
given in Figured 


Theorem 4.4. Consider the GGL hazard function in (JfW. Then, for a > 1, the hazard 
function is increasing and for 0 < a < 1 , is decreasing and bathtub shaped. 

Proof. The proof is omitted, since 0 < 0 < 1 and therefore the proof is similar to the proof of 
Theorem 1. □ 


The hrst and second non-central moments of GGL are 


EiX) = 


afd 




— logfl — 9) ^ 

J n=l fc=0 3=0 


r{k + 1)(7 -I- 7^)2 


-log(l-0)^ F(/i; +1)(7 + 7A:)3 


5 Estimation and inference 

In this section, we will derive the maximum likelihood estimators (MLE) of the unknown 
parameters 0 = {a, (3,^,9)'^ of the GGPS(a,/3, 7 , 0). Also, asymptotic conhdence intervals of 


























these parameters will be derived based on the Fisher information. At the end, we proposed an 
Expectation-Maximization (EM) algorithm for estimating the parameters. 


5.1 MLE for parameters 

Let Xi,...,Xn be an independent random sample, with observed values xi,...,Xn from 
GGPS(q;,/3, 7 , 0) and 0 = (a,/ 3 , 7 , 0 )^ be a parameter vector. The log-likelihood function 
is given by 

n n 

ln = ln{®]x) = nlog( 6 ')-t nlog(Q;/3)-I-n 7 x-F ^log(l - L)-h (a - 1 ) ^log(ti) 

i=l i=l 

n 

+ ^log(C'( 0 tf))-nlog(C( 0 )), 
i=l 

where = 1 — e Therefore, the score function is given by [7(0; x) = (^, 

where 

din 
da 

din 

dp 

din 

d'j 


din 

d9 

The MLE of 0, say 0, is obtained by solving the nonlinear system U{Q;x) = 0. We cannot 
get an explicit form for this nonlinear system of equations and they can be calculated by using 
a numerical method, like the Newton method or the bisection method. 

Eor each element of the power series distributions (geometric, Poisson, logarithmic and 
binomial), we have the following theorems for the MLE of parameters: 

Theorem 5.1. Let gi(a;/?, 7 , 0, x) denote the funetion on RHS of the expression in dSU), 
where p, 7 and 6 are the true values of the parameters. Then, for a given /3 > 0, 7 > 0 and 
9 > 0, the roots 0 /gi(a,/3; 7 , 0, a;) = 0, lies in the interval 

n n 

c'(e) + ^ i=i *=i 





n 

a 


+ '^log{ti) +Y^ 


i=l 


i=l 


9tflog{U)C"{9tf) 
C'{9tf) 


1 ^ 


_ 


n 


n dti n nditf) 

tj. 


2 = 1 


2=1 




2=1 2=1 


Xie 


'yxi^ 


n dti n 

+E ^ 


n 


+ E 


i=i ‘ i=i 

tfC"{9tf) nC'{9) 


C'{9tf 


i=l 


cm 


C{9) 


cm 


(5.1) 

(5.2) 


(5.3) 

(5.4) 











Proof. See Appendix B.l. 


□ 


Theorem 5.2. Let g 2 (/ 3 ; a, 7 , 0, x) denote the function on RHS of the expression in ()5.3p . 
where a, 7 and 9 are the true values of the parameters. Then, the equation g 2 (/ 3 ; a, 7 , 9,x) = 0 
has at least one root. 

Proof. See Appendix B.2. □ 

Theorem 5.3. Aei g 3 ( 0 ; a,/3, 7 , x) denote the function on RHS of the expression in (|5.4p and 
X = n~^ where a, (5 and 7 are the true values of the parameters. 

a) The equation g 3 ( 0 ; a,/3, 7 , x) = 0 has at least one root for all GGG, GGP and GGL distri¬ 
butions if Yl'i=i > i • 

b) If g 3 (p; a, ft, 7 , x) = where p = and p G (0,1) then the equation gz{p', a, /3, 7 , x) = 0 

has at least one root for GGB distribution if X]r=i • 

Proof. See Appendix B.3. □ 


Now, we derive asymptotic confidence intervals for the parameters of GGPS distribution 


It is well-known that under regularity conditions (see 


Casella and Bergen. 


2001 


, Section 10), 


the asymptotic distribution of ^/n(& — 0) is multivariate normal with mean 0 and variance- 
covariance matrix J“^(0), where Jn(0) = lim„_).oo -ln(0)) and 1^(0) is the 4x4 observed 
information matrix, i.e. 
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whose elements are given in Appendix C. Therefore, an 100(1—r/) asymptotic confidence interval 
for each parameter, 0 ^, is given by 

AC R — (0r — ^rj/2 ^ ; 0r 4 “ I, 


where Irr is the (r, r) diagonal element of /„^(0) for r = 1,2,3,4 and Z ^/2 is the quantile ^ 
of the standard normal distribution. 


5.2 EM-algorithm 


The traditional methods to obtain the MLE of parameters are numerical methods by solving the 
equations (|5.in - (l5.4l) . and sensitive to the initial values. Therefore, we develop an Expectation- 
Maximization (EM) algorithm to obtain the MLE of parameters. It is an iterative metho d, 


and is a very powerful tool in handling the incomplete data problem (jPempster et al 


19771). 











We define a hypothetical complete-data distribution with a joint pdf in the form 


5 ((x,2;;0) = S 

where t = 1 — and a,/3, 7 , 0 > 0, x > 0 and z G N. Suppose 

is the current estimate (in the rth iteration) of ©. Then, the E-step of an EM cycle 
requires the expectation of {Z\X;@^'^'^). The pdf of Z given X = x is given by 


g{z\x) 




and since 


c'{e) + ec"{e) = ^ a^ze^-^ + ^ 

2=1 2 = 1 


'^z^azO'' S 

Z = 1 


the expected value of Z\X = x is obtained as 


E{Z\X = x) 




(5.5) 


By using the MLE over 0, with the missing Z’s replaced by their conditional expecta¬ 
tions given above, the M-step of EM cycle is completed. Therefore, the log-likelihood for the 
complete-data is 


n n 

l*n{y,®) OC ^ Zi log(0) -I- n log(a/3) -|- nyx -|- ^ log(l — ti) 

i=\ i=l 

n 

+ '^{zia - l)log(fj) - nlog(C( 6 ')), (5.6) 

i=l 


where y = (x;^;), x = (xi,...,x,i) and 2 ; = (zi,..., Zn)- On differentiation of (|5.6n with 
respect to parameters a, 13, 7 and 6, we obtain the components of the score function, U{y, 0) = 
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da ’ 8/3 ’’ d'Y ’’ 89 ■' ’ 
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From a nonlinear system of equations U{y;@) = 0, we obtain the iterative procedure of the 
EM-algorithm as 
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where ^(i+i) and yb’+i) are found numerically. Here, for i = 1, 2,... , n, we have that 


Za — 1 “b 




where [1 — e 

We can use the results of 


11 ,-s 


Louii 


(jl982l l to obtain the standard errors of the estimators from 
the EM-algorithm. Consider lc{G;x) = E{Ic{&;y)\x), where Ic{Q;y) = is the 

4x4 observed information matrix.If lm{®] x) = Var[U{y, 0)|a;], then, we obtain the observed 
information as 

/(©; x) = lc{&-, x) - lm{®; x). 

The standard errors of the MLEs of the EM-algorithm are the square root of the diagonal 
elements of the/(0; a;). The computation of the se matrices are too long and t edious. Therefore, 


we did not present the details. Reader can see 
these values. 


Mahmoudi and Jafaril (120121 ) how to calculate 


6 Simulation study 

We performed a simulation in order to investigate the proposed estimator of a, /3, 7 and 9 of the 
proposed EM-scheme. We generated 1000 samples of size n from the GGG distribution with 
/3 = 1 and 7 = 0.1. Then, the averages of estimators (AE), standard error of estimators (SEE), 
and averages of standard errors (ASE) of MLEs of the EM-algorithm determined though the 
Fisher information matrix are calculated. The results are given in Table [21 We can find that 






















(i) convergence has been achieved in all cases and this emphasizes the numerical stability of 
the EM-algorithm, 

(ii) the differences between the average estimates and the true values are almost small, 

(hi) the standard errors of the MLEs decrease when the sample size increases. 


Table 2; The average MLEs, standard error of estimators and averages of standard errors for 


the GGG distribution. 



parameter 

AE 

SEE 

ASE 

n 

a 

e 

a 

/? 

7 

e 

a 

/? 

7 

e 

a 

/? 

7 

e 

50 

0.5 

0.2 

0.491 

0.961 

0.149 

0.204 

0.114 

0.338 

0.265 

0.195 

0.173 

0.731 

0.437 

0.782 


0.5 

0.5 

0.540 

0.831 

0.182 

0.389 

0.160 

0.337 

0.260 

0.263 

0.210 

0.689 

0.421 

0.817 


0.5 

0.8 

0.652 

0.735 

0.154 

0.684 

0.304 

0.377 

0.273 

0.335 

0.309 

0.671 

0.422 

0.896 


1.0 

0.2 

0.988 

0.972 

0.129 

0.206 

0.275 

0.319 

0.191 

0.209 

0.356 

0.925 

0.436 

0.939 


1.0 

0.5 

1.027 

0.852 

0.147 

0.402 

0.345 

0.352 

0.226 

0.283 

0.408 

0.873 

0.430 

0.902 


1.0 

0.8 

1.210 

0.711 

0.178 

0.745 

0.553 

0.365 

0.230 

0.342 

0.568 

0.799 

0.433 

0.898 


2.0 

0.2 

1.969 

0.990 

0.084 

0.216 

0.545 

0.305 

0.151 

0.228 

0.766 

1.135 

0.422 

0.902 


2.0 

0.5 

1.957 

0.842 

0.113 

0.487 

0.608 

0.334 

0.192 

0.277 

0.820 

1.061 

0.431 

0.963 


2.0 

0.8 

2.024 

0.713 

0.161 

0.756 

0.715 

0.396 

0.202 

0.353 

1.143 

0.873 

0.402 

0.973 

100 

0.5 

0.2 

0.491 

0.977 

0.081 

0.212 

0.084 

0.252 

0.171 

0.179 

0.125 

0.514 

0.283 

0.561 


0.5 

0.5 

0.528 

0.883 

0.109 

0.549 

0.124 

0.275 

0.178 

0.247 

0.155 

0.504 

0.275 

0.567 


0.5 

0.8 

0.602 

0.793 

0.136 

0.769 

0.215 

0.323 

0.194 

0.299 

0.220 

0.466 

0.259 

0.522 


1.0 

0.2 

0.974 

0.997 

0.102 

0.226 

0.195 

0.242 

0.129 

0.206 

0.251 

0.645 

0.280 

0.767 


1.0 

0.5 

1.030 

0.875 

0.113 

0.517 

0.262 

0.291 

0.155 

0.270 

0.298 

0.651 

0.295 

0.843 


1.0 

0.8 

1.113 

0.899 

0.117 

0.846 

0.412 

0.342 

0.177 

0.331 

0.400 

0.600 

0.287 

0.781 


2.0 

0.2 

1.952 

0.995 

0.138 

0.221 

0.424 

0.237 

0.117 

0.209 

0.524 

0.922 

0.321 

0.992 


2.0 

0.5 

2.004 

0.885 

0.110 

0.518 

0.493 

0.283 

0.131 

0.274 

0.601 

0.873 

0.321 

0.966 


2.0 

0.8 

2.028 

0.981 

0.104 

0.819 

0.605 

0.350 

0.155 

0.339 

0.816 

0.717 

0.289 

0.946 


7 Real examples 

In this Section, we consider two real data sets and fit the Gompertz, GGG, GGP, GGB (with 
m = 5), and GGL distributions. The first data set is negatively skewed, and the second data set 
is positively skewed, and we show that the proposed distributions fit both positively skewed and 
negatively skewed data well. Eor each data, the MLE of parameters (with standard deviations) 
for the distributions are obtained. To test the goodness-of-fit of the distributions, we calculated 
the maximized log-likelihood, the Kolmogorov-Smirnov (K-S) statistic with its respective p- 
value, the AIG (Akaike Information Griterion), AIGG (AIG with correction), BIG (Bayesian 
Information Criterion), CM (Cramer-von Mises statistic) and AD (Anderson-Darling statistic) 
for the six distributions. Here, the significance level is 0.10. To show that the likelihood 
























Table 3: The strengths of glass fibers. 


0.55, 

0.93, 

1.25, 

1.36, 

1.49, 1.52, 

1.58, 

1.61, 1.64, 

1.68, 1.73, 

1.81, 2 . 00 , 

0.74, 

1.04, 1.27, 

1.39, 

1.49, 

1.53, 

1.59, 

1.61, 1 . 66 , 

1 . 68 , 

1.76, 1.82, 

2.01, 0.77, 

1 . 11 , 1.28, 

1.42, 

1.50, 1.54, 

1.60, 

1.62, 

1 . 66 , 

1.69, 

1.76, 1.84, 

2.24, 

0.81, 1.13, 

1.29, 1.48, 

1.50, 1.55, 

1.61, 

1.62, 1 . 66 , 

1.70, 

1.77, 

1.84, 

0.84, 

1.24, 1.30, 

1.48, 

1.51, 1.55, 

1.61, 1.63, 

1.67, 1.70, 

1.78, 

1.89 




Table 4: The phosphorus concentration in 

the leaves. 



0 . 22 , 

0.17, 

0 . 11 , 

0 . 10 , 

0.15, 0.06, 

0.05, 

0.07, 0.12, 

0.09, 0.23, 

0.25, 0.23, 

0.24, 

0 . 20 , 0.08 

0 . 11 , 

0 . 12 , 

0 . 10 , 

0.06, 

0.20, 0.17, 

0 . 20 , 

0 . 11 , 0.16, 

0.09, 0.10, 

0 . 12 , 0 . 12 , 

0 . 10 , 

0.09, 0.17 

0.19, 

0 . 21 , 

0.18, 

0.26, 

0.19, 0.17, 

0.18, 

0.20, 0.24, 

0.19, 0.21, 

0.22, 0.17, 

0.08, 

0.08, 0.06 

0.09, 

0 . 22 , 

0.23, 

0 . 22 , 

0.19, 0.27, 

0.16, 

0.28, 0 . 11 , 

0 . 10 , 0 . 20 , 

0.12, 0.15, 

0.08, 

0.12, 0.09 

0.14, 

0.07, 

0.09, 

0.05, 

0.06, 0 . 11 , 

0.16, 

0.20, 0.25, 

0.16, 0.13, 

0 . 11 , 0 . 11 , 

0 . 11 , 

0.08, 0.22 

0 . 11 , 

0.13, 

0 . 12 , 

0.15, 

0 . 12 , 0 . 11 , 

0 . 11 , 

0.15, 0.10, 

0.15, 0.17, 

0.14, 0.12, 

0.18, 

0.14, 0.18 

0.13, 

0 . 12 , 

0.14, 

0.09, 

0.10, 0.13, 

0.09, 

0 . 11 , 0 . 11 , 

0.14, 0.07, 

0.07, 0.19, 

0.17, 

0.18, 0.16 

0.19, 

0.15, 

0.07, 

0.09, 

0.17, 0.10, 

0.08, 

0.15, 0.21, 

0.16, 0.08, 

0 . 10 , 0.06, 

0.08, 

0.12, 0.13 


equations have a unique solution in the parameters, we plot the profile log-likelihood functions 
of /3, 7 , a and Q for the six distributions. 

First, we con s ider t he data consisting of the strengths of 1.5 cm glass fibers given in 


Smith and Navlon (llQSTIi and measured at t 


data is also studied by 


le b 


Barreto-Souza et al. 


ational Physical Laboratory, England. This 


( 201C ) and is given in Table [3l 


The results are given in Table [5] and show that the GGG distribution yields the best fit 
among the GGP, GGB, GGL, GG and Gompertz distributions. Also, the GGG, GGP, and 
GGB distribution are better than GG distribution. The plots of the pdfs (together with the 
data histogram) and cdfs in Figure [7] confirm this conclusion. Figures [9] show the profile log- 

likelihood functions of /3, 7 , a and d for the six distrib utions. _ 

As a second example, we consider a data set from iFonseca and Francal (j2007l i. who stud¬ 


ied the soil fertility influence and the characterization of the biologic fixation of N 2 for the 
Dimorphandra wilsonii rizz growth. For 128 plants, they made meas u res o f the phosphorus 
concentration in the leaves. This data is also studied by 


Silva et al. 


(l2013l i and is given in 


Table 01 Figures [TO] show the profile log-likelihood functions of (3, 7 , a and 9 for the six 
distributions. 

The results are given in Table [ 6 l Since the estimation of parameter 9 for GGP, GGB, 
and GGL is close to zero, the estimations of parameters for these distributions are equal to 
the estimations of parameters for GG distribution. In fact. The limiting distribution of GGPS 
when 0 —)• 0^ is a GG distribution (see Proposition [2|) . Therefore, the value of maximized log- 
likelihood, log(L), are equal for these four distributions. The plots of the pdfs (together with 



























the data histogram) and cdfs in Figure [8] confirm these conclusions. Note that the estimations 
of parameters for GGG distribution are not equal to the estimations of parameters for GG 
distribution. But the log(L)’s are equal for these distributions. However, from Table [6] also we 
can conclude that the GG distribution is simpler than other distribution because it has three 
parameter but GGG, GGP, GGB, and GGL have four parameter. Note that GG is a special 
case of GGPS family. 
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Figure 7: Plots (pdf and cdf) of fitted Gompertz, generalized Gompertz, GGG, GGP, GGB 
and GGL distributions for the first data set. 
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Figure 8: Plots (pdf and cdf) of fitted Gompertz, generalized Gompertz, GGG, GGP, GGB 
and GGL distributions for the second data set. 






























Table 5; Parameter estimates (with std.), K-S statistic, p-value, AIC, AICC and BIC for the 
hrst data set,_ 


Distribution 

Gompertz 

GG 

GGG 

GGP 

GGB 

GGL 

/3 

0.0088 

0.0356 

0.7320 

0.1404 

0.1032 

0.1705 

s.e.(/3) 

0.0043 

0.0402 

0.2484 

0.1368 

0.1039 

0.2571 

7 

3.6474 

2.8834 

1.3499 

2.1928 

2.3489 

2.1502 

s.e.{j) 

0.2992 

0.6346 

0.3290 

0.5867 

0.6010 

0.7667 

a 

— 

1.6059 

2.1853 

1.6205 

1.5999 

2.2177 

s.e.{a) 

— 

0.6540 

1.2470 

0.9998 

0.9081 

1.3905 

e 

— 

— 

0.9546 

2.6078 

0.6558 

0.8890 

s.e.{6) 

— 

— 

0.0556 

1.6313 

0.5689 

0.2467 

-log(L) 

14.8081 

14.1452 

12.0529 

13.0486 

13.2670 

13.6398 

K-S 

0.1268 

0.1318 

0.0993 

0.1131 

0.1167 

0.1353 

p-value 

0.2636 

0.2239 

0.5629 

0.3961 

0.3570 

0.1992 

AIC 

33.6162 

34.2904 

32.1059 

34.0971 

34.5340 

35.2796 

AICC 

33.8162 

34.6972 

32.7956 

34.78678 

35.2236 

35.9692 

BIC 

37.9025 

40.7198 

40.6784 

42.6696 

43.1065 

43.8521 

CM 

0.1616 

0.1564 

0.0792 

0.1088 

0.1172 

0.1542 

AD 

0.9062 

0.8864 

0.5103 

0.6605 

0.7012 

0.8331 


Table 6; Parameter estimates (with std.), K-S statistic, p-value, AIC, AICC and BIC for the 
second data set._ 


Distribution 

Gompertz 

GG 

GGG 

GGP 

GGB 

GGL 

/3 

1.3231 

13.3618 

10.8956 

13.3618 

13.3618 

13.3618 

s.e.0) 

0.2797 

4.5733 

8.4255 

5.8585 

6.3389 

7.3125 

7 


3.1500 

4.0158 

3.1500 

3.1500 

3.1500 

s.e.{j) 


2.1865 

3.6448 

2.4884 

2.6095 

2.5024 

a 

— 

6.0906 

5.4236 

6.0906 

6.0906 

6.0905 

s.e.(d) 

— 

2.4312 

2.8804 

2.6246 

2.7055 

2.8251 

e 

— 

— 

-0.3429 

1.0 X 10"® 

1.0 X 10"® 

1.0 X 10"® 

s.e.{6) 

— 

— 

1.2797 

0.8151 

0.2441 

0.6333 

- log(T) 

-184.5971 

-197.1326 

-197.1811 

-197.1326 

-197.1326 

-197.1326 

K-S 

0.1169 

0.0923 

0.0898 

0.0923 

0.0923 

0.0923 

p-value 

0.06022 

0.2259 

0.2523 

0.2259 

0.2259 

0.2259 

AIG 

-365.1943 

-388.2653 

-386.3623 

-386.2653 

-386.2653 

-386.2653 

AIGG 

-365.0983 

-388.0717 

-386.0371 

-385.9401 

-385.9401 

-385.9401 

BIG 

-359.4902 

-379.7092 

-374.9542 

-374.8571 

-374.8571 

-374.8571 

CM 

0.3343 

0.1379 

0.1356 

0.1379 

0.1379 

0.1379 

AD 

2.3291 

0.7730 

0.7646 

0.7730 

0.7730 

0.7730 
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Figure 9: The profile log-likelihood functions for Gompertz, generalized Gompertz, GGG, GGP, 
GGB and GGL distributions for the first data set. 
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Figure 10: The profile log-likelihood functions for Gompertz, generalized Gompertz, GGG, 


GGP, GGB and GGL distributions for the second data set. 

















Appendix 


We demonstrate those parameter intervals for which the hazard function is decreasing, i ncreas¬ 
i ng an d bathtub shaped, and in order to do so, we follow closely a theorem given by 


Glaser 


(j 19801 1. Define the function t{x) = where f'{x) denotes the first derivative of f{x) in 


(12.3p . To simplify, we consider u = 1 — exp(^(e'>'^ “ !))• 


A.l 

Consider the GGG hazard function in (14.2h . then we define 

/ \ 1 — a 2a9u°^~^ 

"W = 7(T = — + 

If a > 1, then t '{ u ) > 0, and h{.) is an increasing function. If 0 < a < 1, then 

, , 2a6'^ 1 

lim r (u) = -oo, lim r (u) = + (a - 1)(1 - > 0. 

Since the limits have different signs, the equation t '{ u ) = 0 has at least one root. Also, we 
can show that t "{ u ) > 0. Therefore, the equation t '{ u ) = 0 has one root. Thus the hazard 
function is decreasing and bathtub shaped in this case. 


A.2 


The GGP hazard rate is given by h{u) = 6a(5u'^ — e®““). We define = log[/i(n)]. 

Then, its hrst derivative is 


/ / X a — I ^ n M 

Tj {u) =-h aOe^ 


OL—1 


U 


g6» _ 


It is clearly for a > 1, r]'{u) > 0 and h{u) is increasing function. If 0 < a < 1, then 


lim rj'{u) = —oo, lim rj'{u) = 0, 

So the equation t'{u) = 0 has at least one root. Also, we can show that r"(u) > 0. It implies 
that equation r]'(u) = 0 has a one root and the hazard rate increase and bathtub shaped. 


B. 


B.l 


Let to,(a) = Eti = W Eti log(C'(«i“)). For GGG, 


wi{a) = 26 


2=1 


tf log tj dwi{a) 

1 — 9tf ’ da 


» Li _ g^al 


2=1 
















For GGP, 


For GGL, 


rt;i(a) = O'^tflogU, 


i=l 


wi{a) = 9 


dwi{a) 
da 

tf log ti dwi{a) 


> 0 . 


2=1 




2=1 


1 - 6»i“ ’ 


da 


2=1 


1 - 9t^ 


For GGB, 


u,(a) = (m - 1)9 V = {.» - 1)9 E ‘f [^1" > 0. 


i=l 


2=1 


Therefore, rt;i(a) is strictly increasing in a and 


lim gi{a;/3,-f,9,x) = oo, lim gi(a;/?, 7 ,0, x) = y^logiU). 

Q->- 0 + a^oo 

1=1 

Also, 

n ^ ^ ji 6C"(6) ^ ^ 

gi(a;^,7,6',x) < - + ^log{ti), gi{a;/ 3 ,j,e,x) > - + ( + l)^log{ti). 

“ i=i a o 

Hence, gi(a; 7 , 0, x) < 0 when ^ + 1]?=! log(^*) < 0i and gi(a;/3, 7 , 6 », x) > 0 when | + 
Sr=i log(^*) > 0 - The proof is completed. 

B.2 

It can be easily shown that 

-1 ” 

lim g 2 (/ 3 ;a, 7 , 6 »,x) = 00 , lim g 2 (/ 3 ; a, 7 , 0 , x) = —- 1). 

/3-)^0+ /9-)^oo 7 ^ 

^=l 

Since the limits have different signs, the equation g 2 (/ 3 ; a, 7 ,0, x) = 0 has at least one root with 
respect to /3 for fixed values a, 7 and 9. The proof is completed. 


B.3 

a) For GGP, it is clear that 


limg 3 ( 6 »;a,/ 3 , 7 ,x) = lim g 3 ( 6 »; a,/?, 7 , x) = - 00 . 

e-l-O ^^ I d^oo 

i=l 

Therefore, the equation g 3 ( 0 ; a, /3, 7 , x) = 0 has at least one root for 0 > 0, if Yll=i “ f > 0 

or > §• 

b) For GGG, it is clear that 


lim g 3 ( 6 »;a,/ 3 , 7 ,x) = - 00 , 
0^00 


n 


lim g 3 (( 9 ; 0 ,^, 7 , x) 
0-)-O+ 


= -n + 2Y,tt 

2 = 1 



Therefore, the equation g3{9, (3, 7 , x) = 0 has at least one root for 0 < 0 < 1, if —n+2 > 

0 or > §• 

For GGL, it is clear that 


liing 3 ( 6 »;a,/ 3 , 7 ,x) = Vtf - lim g 3 ( 6 '; a,/3, 7 , x) = - 00 . 


2=1 


2’ 0^1 ‘ 


Therefore, the equation g 3 ( 0 ; a, (3, 7 , x) = 0 has at least one root for 0 < 0 < 1, if ^11=1 “ § > 

Oor >§• 

For GGB, it is clear that 


liing 3 (p;a,/ 3 , 7 ,x) = 1) - , lini g 3 (p; a, 7 , x) ^ ^ 1 

p—>-0 ^^ z p —>0 ^ 

i=l i=l 




Therefore, the equation g 3 (p; a, /3, 7 , x) = 0 has at least one root for 0 < p < 1, if XlILi tfim — 

m- 

~T 


_ n(m-l) 


> 0 and E?.i < 0 or E?.i <f > ? and E?=i *r" > iSk- 


c. 

Consider t* = 1 — e 7 ^*^^ ' Then, the elements of 4 x 4 observed information matrix In{Q) 
are given by 


laa or 9 — 

oa‘^ a 


n/^ < ^ I 

2=1 






Iaj3 — 


d'^L 


dad 13 ^ ■ 7 


"■ g7d;i _ 2 0 ”■ 


E 


+ -111(0 iog(t.) +1)§^ 

^ 2=1 ^ ^ 


^ 0:7 — 


+Q;6»tf log(ti) 

d^ln 
dadj 


C'"{9tf)C'{etf) - {C''{9tf)f 




{C'{9tf)f 

e7dir(7x. _i) + i ep 


2=1 


7" 


] + 4 + l][(alog(t,) + 1 ) 


3-a9 


I(3f} = 


Ip'y — 


+a9ti log(ti) 

dX 

dad9 


C"'{9tf)C'{9tf) - {C''{9tf)) 


^ ^=l 

ll(Q+a\\2 


C''{9tt) 

C'{9tf) 




2=1 


ic'rnw 

C"{9tf) C'''{9tf)C’{9tf) - {C"{9t<t)f 

C'{9tf) * 


{C'{9tf)r 


—n 


dX _ 

dl3^ /32 

dX {a - 2 ) 


+ 9a‘^'^tj 


^^^a^tf)C'{9tf) - iC"{9tf)y 


2=1 

n 


7 ^ X'{9tf) 


{C'{9tf)Y 


dpd'y 


E(*""-(7a^. -D + D + aOY. - ') + ‘>1?^ 


2=1 

^ C'''i9tf)C'i9X-iC''i9tf)y 

X ^ X{9tf)W 



Ifie = 


-^77 ~ 


lO'y 

lee 


dHn 

d/3de 


E*: 

2 = 1 


l(aj.a\'\2 


ic'ietf)) 


rfil OR JE —2 q.'^p'Y^i 

- 1) + 1] + (a - - 1) + 1) + ^1 


2=1 


-2 


2=1 


■ 


/-yO 




7 metf) 




7 ' 


+ ^(e^"*(7x.-l) + l) 
T 

dHn 
d9dj 


C"(0tf) 


( C "( 02 f )))2 


_ 1) +1,1^ + - (a”(«.“))^ 

i=l ^ 


C"(0tf 


502 02 


-n ^.2a. C'"{9tf)C'{0tf)-{C"{0t<l)r 
ft2 Z-,- (riltQ+a\\2 J 

2=1 


(C'(0tf))2 


(C'(02f))2 

C''{0)C'{0) - (C'(0))2 
(C"(0))2 ^ 
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